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The evolution of a two-dimensional driven lattice-gas model is studied on an L x x L y lattice. 
Scaling arguments and extensive numerical simulations are used to show that starting from random 
initial configuration the model evolves via two stages: (a) an early stage in which alternating stripes 
of particles and vacancies are formed along the direction y of the driving field, and (b) a stripe 
coarsening stage, in which the number of stripes is reduced and their average width increases. The 
number of stripes formed at the end of the first stage is shown to be a function of L x /Ly, with 
<f> ~ 0.2. Thus, depending on this parameter, the resulting state could be either single or multi 
striped. In the second, stripe coarsening stage, the coarsening time is found to be proportional 
to L y , becoming infinitely long in the thermodynamic limit. This implies that the multi striped 
state is thermodynamically stable. The results put previous studies of the model in a more general 
framework. 
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I. INTRODUCTION 

Driven diffusive systems have been extensively studied 
in recent years. They serve as a fruitful framework for 
studying the statistical mechanics of systems far from 
thermal equilibrium. Driven by an external field these 
systems reach a steady state with a non-vanishing cur- 
rent and as such do not satisfy detailed balance. Studies 
of these models have revealed many differences between 
such systems and systems in thermal equilibrium. For ex- 
ample, several one-dimensional driven diffusive systems 
with local dynamics exhibit long range order and spon- 
taneous symmetry breaking. Such phenomena can not 
occur in thermal equilibrium when the interactions are 
short ranged fl]||. 

Many studies of driven diffusive systems have focused 
on a driven lattice gas (Ising) model. The model was in- 
troduced by Katz, Lebowitz, and Spohn ||] and is often 
referred to as the "standard model" . In d = 2 dimen- 
sions the model is defined on an L x x L y lattice. Each 
of the lattice sites i, is either occupied by a particle or is 
vacant. A macroscopic configuration is characterized by 
a set of occupation numbers {rii} where ni = 0, 1 repre- 
sents a vacant or an occupied site, respectively. Usually 
the model is studied with an equal number of occupied 
and vacant sites. An energy TL = — n i n j is associ- 
ated with each configuration. Here the sum is over (ij) 
nearest neighbor sites. The energy represents an attrac- 
tive interaction between the particles. An external drive 
is introduced through a field E which biases the motion 
of the particles in the —y direction. Imposing periodic 
boundary conditions in this direction results in a current 
of particles through the system along the field direction. 
Specifically the dynamics of the model is defined through 
the exchange of nearest neighbor particles with a rate 

W = mm{l,exp(-0AH-EAy)}. (1) 
Here (3 is an inverse temperature-like parameter, and 



Ay = (—1,0, 1) for a particle attempting to hop along, 
orthogonal to, or against the direction of the driving field. 
The energy difference between the two configurations af- 
ter and before the particle exchange is denoted by AH. 

The model has been studied extensively for nearly two 
decades [Q. Monte-Carlo simulations suggest that the 
(T, E) phase diagram of the model is composed of two 
phases: A high temperature disordered phase in which 
the particle density is homogeneous, and a low temper- 
ature phase in which the system orders and phase sepa- 
rates into high density and low density regimes. It was 
found that in this phase the particles evolve towards a 
striped structure parallel to the direction of the driving 
field. Numerical studies indicate that a slow coarsening 
takes place in this state As the magnitude of the 

driving field is increased, the transition temperature be- 
tween the two phases increases and saturates at about 
1.41 To ||, where To is the Onsager temperature corre- 
sponding to E = 0. 

Recent Monte-Carlo simulations of this model suggest 
that the evolution of the striped phase is rather complex. 
For a square system the stripes are found not to coarsen 
in the thermodynamic limit, yielding a multi-striped or- 
dered state. This phase was termed extraordinary or 
"stringy" ju^. On the other hand systems with large 
aspect ration, L y 3> L x , were found to evolve toward a 
single stripe phase. 

In order to get a better understanding of the nature 
of the ordered phase of the driven lattice gas model we 
carry out in this paper a finite size scaling analysis of the 
evolution process starting from a fully disordered state. 
We find that the model evolves via two stages: 

(a) an early stripe formation stage in which stripes are 
formed from the initially disordered state; and 

(b) a stripe coarsening stage in which the multi-stripe 
configuration formed in the early stage coarsen by reduc- 
ing the number of stripes and increasing their average 
width. A typical evolution of such a system is shown in 
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FIG. 1. A typical evolution of a system of size 
L x = L y — 100 from a random initial condition. Configura- 
tions from times (a) 50, (b) 2000 and (c) 500,000 Monte-Carlo 
sweeps are shown. Here (5 — 2 and E = oo. One can clearly 
observe the two steps of the coarsening process described in 
the text. 



Fig. 

Our studies yield two main results: 

(1) The number of stripes which are formed at the end of 
the initial stripe formation stage strongly depends on the 
aspect ratio of the system. In particular we find that the 
number of stripes m scales asm~ L x /L^, with <f> ~ 0.2. 
This implies that for narrow systems (L x /L^ < 1) a sin- 
gle stripe is formed at the end of the first stage, while 
for wide systems (L x /L^ ^> 1) the resulting structure is 
multi-striped. 

(2) Simple arguments are presented to show that dur- 
ing the stripe coarsening stage the average width of the 
stripes grows with time as (t/Ly) 1 / 3 . This behavior 
is verified by extensive numerical simulations. There- 
fore, the coarsening of the stripes becomes slower as the 
system size in the direction of the drive is increased. 
This implies that in the thermodynamic limit a multiple 
striped configuration is in fact stable. We note that sim- 
ilar phenomena of arrested striped configurations have 
been observed in previous studies of coarsening of other 
models with striped structures perpendicular to the di- 
rection of the drive Jll]] . 

The paper is organized as follows: In Section II the 
stripe formation stage is discussed. Section III considers 
the stripe coarsening stage. We end with a summary 
and discussion of the implications of our results to other 
related works in Section IV. 



II. THE STRIPE FORMATION STAGE 

The evolution of the driven lattice gas model in the 
early stripe formation stage has received some atten- 
tion |6| Numerical simulations indicate that the do- 
main growth process which takes place in this stage is 
highly anisotropic. The typical domain size in the di- 
rection of the drive and the direction perpendicular to it 
grow differently. In particular it has been observed [^]]7]] 
that the typical domain size parallel to the drive grows 
roughly as l y ~ i^u with tpn ~ 1, while the typical do- 



main size perpendicular to the drive grows roughly as 
~ t v± , cpx — 0.2. This behavior is very different 
from that of a non-driven system evolving towards equi- 
librium. It is well known that in such a system, when the 
dynamics is conserving, as is the case here, the average 
linear domain size £ grows as t 1 / 3 U). The difference 
in behavior is due to the inherent anisotropy induced by 
the drive. 

The number of stripes formed in the system at the 
end of the stripe formation stage can be estimated using 
the results described above. For a stripe to form in the 
system the size of a domain along the direction of the 
drive t y (t) must be of the order of the system size L y . 
Since l y (t) ~ the time for this to occur t s scales as 
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At this time the typical domain size perpen- 



dicular to the drive is 



Thus, the number of stripes formed, m, scales as 
L x L x 

~ if 



m 



(2) 



(3) 



where <f> = ip±/ip^. Using the estimates for the exponents 
ip±_ and ip« one has <f> ~ 0.2. 

Specifically, for a square system, where L x = L y = L, 
Eq. H implies m ~ L 1 ^^. Since <ft < I we find that the 
number of stripes grows as the system size is increased, 
and one always reaches a multi-striped state. The stripe 
density, m/L x , vanishes in the thermodynamic limit. In 
fact, this result is independent of the exact value of the 
exponent <p as long as tpu > tpj_. 

To verify these results Monte-Carlo simulations are 
performed for various system sizes, starting from a ran- 
dom initial condition. The Monte-Carlo procedure we 
use is standard: At each time step a pair of neighbor- 
ing sites is chosen randomly and updated according to 
the rate W given in (|l|). Throughout the paper we use 
E = oo and (3 — 2, for which the system is ordered. 
We have checked that the main features of this study are 
unchanged for other values of the parameters as long as 
the system is in the ordered phase. We first verify the 
growth law of t s with L y . In order to evaluate t s the 
equal-time correlation of two sites at a distance L y /2 in 
the drive direction is measured and averaged over the 
sample. The time t s is estimated by the time at which 
the average measured correlation reaches the value of 0.4. 
The results are shown in Fig. ||(a) . One can see that the 
behavior of t s with L y is consistent with tpa ~ 1. 

The number of stripes initially formed in the system 
is estimated by performing a Fourier transform of the 
density in the x direction and locating its first peak at 
a non-zero wave length. This procedure is repeated 40 
times for each system size. For simplicity we consider 
only square systems. In Fig. we plot the location 
of the peak for square systems as a function of L. The 
fitted exponent for slightly over a decade of system sizes 
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FIG. 2. (a) The stripe formation time t s (in Monte-Carlo 
sweeps) plotted against L y , for systems with L x = 100. The 
behavior is consistent with <pu ~ 1. (b) The number of stripes 
formed at the end of the first stage in square systems of var- 
ious sizes L = L x — L y , plotted on a log-log scale. The 
straight line corresponds tom~ L 0,82 . 



gives 4> ~ 0.18 ± 0.05 which fits rather well with the 
values predicted by the argument <fi w 0.2. The linear 
dependence of m on L x for non square systems is also 
verified through simulations which are not shown in this 
paper. 

In general, the number of stripes also depends on the 
magnitude of the driving field and the temperature. One 
can write, based on Eq. |[ m — AL X /L^, where the am- 
plitude A(E, f3) is introduced. We find that the ampli- 
tude is an increasing function of both the magnitude of 
the driving field E and the inverse temperature (3. We 
note that when > AL X the number of stripes already 
at the end of the first stage is expected to be one. 



III. THE STRIPE COARSENING STAGE 

We now turn to the second stage of the ordering pro- 
cess, namely the coarsening of the stripes formed in the 
early stage. We present a simple argument suggesting 
that the time t in which stripes coarsen scales linearly 
with the system size parallel to the drive, L y . More- 
over, we show that the average stripe width £(t) scales 
with time as (t/Ly) 1 ^ 3 . The fact that the characteristic 
time associated with the coarsening process scales lin- 
early with L y implies that in the thermodynamic limit, 
where L y — > oo, the coarsening time becomes infinite 
and thus the multi-stripe structure exists as a thermo- 
dynamically stable state. Our argument relies on two 



main features of the driven system: (i) the fact that the 
ordered domains, namely the stripes, are of the size of 
the system, and (ii) the smoothness of the domain walls 
bounding the stripes. This last feature has been shown 
to be a result of the drive fla-llq]. In contrast to the 
non-driven two-dimensional Ising model, where the do- 
main walls are rough, here the driving field makes the 
domain walls smooth. 

We proceed by considering a striped state composed of 
alternating stripes of particles and vacancies with aver- 
age width £. Neighboring stripes of particles interact with 
each other by an exchange of particles. Since the bound- 
aries of the stripes are smooth, the lateral distance that 
particles have to travel in order to move from one stripe 
to the other is of the order of £. To estimate the coarsen- 
ing time we assume that within a stripe of vacancies the 
density of particles is low enough so that the particles 
may be considered as non-interacting. This assumption 
is qualitatively supported by the configurations observed 
in simulations (see, e.g, Fig. [I]). When a particle reaches 
the boundary it is absorbed in the neighboring particle 
stripe. Thus the lateral motion of the particles within a 
stripe of vacancies can be considered as a one-dimensional 
random walk in the x direction with two absorbing walls 
located at x = and x = £. This problem is known as 
the gambler's ruin problem JL7| . The probability of such 
a particle to move from to £ is given by p(£) ~ l/l. 

For the width of a stripe to decrease by one lattice 
spacing it has to lose L y particles. Due to the right-left 
symmetry of the problem, the particle currents from one 
stripe to the other are balanced on average. Therefore a 
net transfer of particles from one stripe to another is only 
due to fluctuations in the lateral current. The net excess 
in the number of particles transferred at a time interval 
t is then proportional to J L y p(£)t. For one stripe to 
shrink and disappear £L y particles must be transferred 
so that 



^L yP (£)t~£L y 



(4) 



Combining this result with p{£) ~ l/£ one finds that the 
average stripe width in the system grows as 
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(5) 



This suggests that the coarsening time scales with L y , 
yielding a stable striped structure in the thermodynamic 
limit. 

The scaling form (^|) may be verified numerically by 
studying the two point particle-particle correlation func- 
tion. To carry out this analysis we note that in an 
isotropic system without a driving field, the coarsen- 
ing process is characterized by a single length scale £(i), 
which could be the linear size of the growing domains. In 
this case the two point particle-particle correlation func- 
tion obeys a scaling form p2[ 
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where r is the distance between two points. Driven sys- 
tems, on the other hand, are non-isotropic, and corre- 
lations along the drive and perpendicular to it behave 
differently. The typical length scale perpendicular to the 
drive is given by (g). Thus we expect the correlation 
function in the x direction to be of the form 



C±(x,t) = g x 



The asymptotic behavior of g±(z) for z — > is expected 
to obey Porod's law, which states that g±(z) = 1/2 — 
r\z with some constant r\. For z — > 00 one should have 
ffxCO-1/4. 

We now turn to describe numerical studies which sup- 
port our results. Note that the scaling variable in (^) in- 
volves three parameters. All three parameters are varied 
in our numerical studies. This is a demanding computa- 
tional task, and a good collapse of the data is a strong 
conformation of the scaling analysis. In these studies an 
initial striped configuration along the drive direction y is 
considered, and its evolution is simulated. The two point 
correlation function Cj_(x,t) is then calculated, and is 
shown to obey the scaling form (|7|). The widths of the 
stripes in the initial configuration are randomly chosen 
from a Poisson distribution with a mean width £q. The 
simulations are performed for lattices of three different 
sizes: 960 x 8, 800 x 16 and 960 x 32. We consider several 
values of L x to demonstrate that this parameter does not 
play an important role in the process. The mean width 
of the stripes in the initial configurations is taken to be 
£q = 4. The two-point particle-particle correlation func- 
tion C±(x,t) is measured and averaged over 110, 75 and 



54 simulations for L v 



16 and 32, respectively. 



The scaling form (Q) suggests that data collapse should 
take place with respect to the three variables x, t and L y . 
This collapse is checked in two steps. First we consider 
L y = 8 and show that C_l is a function of x/t 1 ! 3 as 
expected. This is demonstrated in the inset of Fig. ||. 
Similar results are obtained for the other system sizes as 
well. 

Next, we verify the full scaling form (0). In Fig. || 
correlation functions for the three different system sizes 
are plotted. For each system size the correlation function 
is evaluated for arbitrarily chosen t and the data is then 
plotted as a function of the scaling variable xj (t/Ly) 1 ! 3 . 
Again, the quality of the data collapse supports our main 
result. Although computation time limits us to a rela- 
tively small systems, we believe the quality of the data 
backs our scaling argument. 
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FIG. 3. The two-point particle-particle correlation function 



C± is plotted as a function of the scaling variable x/(t/L l 



)V 3 ; 

for systems of size 960 x 8 (marked by x), 800 x 16 (o), 
960 x 32 (+). The times of measurement are chosen arbi- 
trarily (t = 1,3,8 x 10 6 Monte-Carlo sweeps, respectively). 
In the inset, C± as a function of x/t 1 / 3 is shown for a sys- 
tem of size 960 x 8 and times t = 0.2(x), 1(a), 2(+) xlO 6 
Monte-Carlo sweeps. 



IV. DISCUSSION 

The evolution of the driven lattice gas model was con- 
sidered starting from a random initial configuration. We 
have shown, using simple scaling arguments and exten- 
sive numerical simulations, that the evolution proceeds 
via two stages: an early, stripe formation stage in which 
stripes of the size of the system are formed, followed by 
a second stage in which the stripes coarsen. While the 
first stage lasts t s ~ L y , the system evolves towards a 
single stripe configuration in the second stage at a time 
of order ~ L 3 L y . This is a result of the fact that the 
typical width of stripes in the coarsening stage scales 
with time as £(t) ~ (t/Ly) 1 / 3 . This result indicates that 
the coarsening time of multi-striped configurations scales 
with the system length L V) suggesting that these configu- 
rations exist as stable states in the thermodynamic limit 
L y — ► 00. 

Thus, starting from a random initial configuration, the 
system evolves to one of two types of states, depending 
on its aspect ratio. For L x /L^ < 1 (<j) ~ 0.2) the stripe 
formation stage leads directly to a single stripe state, 
while for L x /L^ ^> 1 multi-striped states are reached. 
The coarsening process of these states proceeds with a 
time scale proportional to L y . 

These results put in a more general framework pre- 
vious studies of this model which considered either the 
early stages [3-|| of the evolution or the nature of the 
steady state \4}. A recent study of a square system has 
shown |1Q| that a multi striped state (termed " stringy" ) is 
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reached from a random initial condition. It was suggested 
that this state is stable. Our studies indicate that this is 
indeed the case for an infinitely large system. However we 
expect a finite system to coarsen to a single striped state 
at a time of the order of L\L y . The fact that the steady 
state of a system with a small aspect ration, L x /L y , was 
found to be composed of a single stripe is consistent with 
our scaling picture. 

Finally, note that the slow coarsening of the stripes is a 
direct consequence of stripes spanning the entire system. 
This is a result of the existence of the drive, and is ex- 
pected to be valid also in higher dimensions. It would be 
interesting to study such processes in high-dimensional 
systems. We note, however, that already in two dimen- 
sions the computational effort was considerable. 
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